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Abstract 

We develop a general algorithm for describing angular decay distributions in 
cascade decay chains of arbitrary length. The general algorithm is used to study 
joint angular decay distributions for the cascade decay B — » > > Dir) + 

7r) + W(^> Iv) where the D** is a generic P-wave charm meson state. Lepton mass 
effects are fully incorporated The joint angular decay distribution depends on 43 
independently measurable decay parameters if the spin parity of the D** is 1 + and 
on 48 decay parameters if the spin parity of the D** is 2 + . We give expressions 
for these decay parameters in terms of the helicity amplitudes of the two-body 
decay processes. An absolute prediction for all the parameters is presented in the 
framework of the heavy quark effective theory. A method for obtaining the helicity 
amplitudes from measured joint angular distributions is suggested. 
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1 Introduction 



Semileptonic decays of B mesons have been much studied both experimentally and the- 
oretically. For recent reviews cf. e.g. j l| for experiment and @ for theory. On the 
experimental side ARGUS [|J and CLEO f|, ^| find that the ground state to ground state 
decays[] B — > Dlv and B — > D*lv do not saturate the total semileptonic decay width. For 
the contribution of the ground state to excited state transitions ARGUS finds (40 ± 10) 
per cent fU and CLEO has (32 ±5) per cent ||. Not much is known experimentally about 
the nature of these other channels. It is, however, natural to expect that the P-wave reso- 
nances D** should play an important role in them. The Particle Data Group || lists two 
P-wave candidates, the 2 + resonance ^2(2460) and the 1 + resonance Pi (2420). They can 
be plausibly interpreted as composed of a c quark almost at rest in the meson rest frame 
and a P-wave light antiquark with total angular momentum j = 3/2 0. Moreover, one 
expects two more D** resonances (0 + , 1 + ) with the P wave antiquark in the j = 1/2 
state. 

The theory of the reactions B — > D**(— > . . .) + W(— > Iv) and of the subsequent 
D** decays has been discussed by a number of authors [7-14]. From angular momentum 
and parity conservation the allowed D** decays are P**(0 + ) — > (D + ir) s , P**(l + ) — > 
(D*n)sorD, D**(2+) -> (D*n) D , and P**(2+) -> (Dn) D , where the subscript denotes the 
partial wave of the decay channel. In the limit m c — > oo one finds that the P**(l + , j = |) 
decays only into (D*n) D and the P**(l + , j = ~) only into (D*n) s 0. This follows from 
angular momentum conservation in the decay of the P-wave antiquark into an 5*-wave 
antiquark and a pion. Moreover, for the P**(2 + ) decays the ratio of probabilities for the 
channels D*n and Pvr is 3:2 in this limit J7|, g§ 

Testing how well these predictions apply to mesons containing finite mass c-quarks is 
very interesting, but difficult for a variety of reasons. In particular, on the theoretical 
side, one has to introduce somewhat arbitrary kinematical corrections (cf. e.g. 0) when 
comparing different decays. The heavy quark symmetry predictions for the decay distri- 
butions in a given channel are more reliable than a comparison of different decays. For 
example, the predictions of the nonrelativistic quark model have been amazingly success- 
ful when they were applied to a given decay channel. However, the angular distribution 
for the decays D** — ► D*n does not contain enough information by itself. Therefore, e.g. 
Ming-Lu, Wise and Isgur |7|] proposed to study in addition to standard decay distributions 
also the decays into longitudinally polarized P*'s. 

We generalize this approach to a complete analysis of the angular decay distributions 
for the full cascade decay 

B -> P**(-> D*(^ Dti) + tt) +W(-> Iv) (1) 

These angular decay distributions depend on three pairs of spherical angles defined by: (i) 
the D* in the D** rest frame, (ii) the D in the D* rest frame and (iii) the lepton in the rest 
frame of the virtual iy-boson. Other popular distributions, like the energy distributions 
for the lepton in the B rest frame or angular decay distributions for the D* in special spin 
states can be expressed by these basic distributions (for given masses of the intermediate 
particles^]). The angular distribution of the D** in the B rest frame is of no interest, since 

^^Here and in the following when bars over the antiparticles are omitted, we mean the pair of charge 
conjugated decays 

2 Here and in the following particle stands for a particle, a resonance or the virtual VF-boson. 
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it has to be spherically symmetric by angular momentum conservation. Lepton mass 
effects are fully incorporated in our approach. Thus our results are applicable also for the 
decays involving the r-lepton. 

The angular decay distributions contain many independently measurable parameters: 
43 for each of the two D**(l + ) cases and 48 for the D**(2 + ) case. All these parameters 
can be expressed in a model independent way by the helicity amplitudes. The overall 
normalization and the overall phase of the helicity amplitudes drop out and this reduces 
the number of real parameters to eight in each of the 1 + cases and to six in the 2 + 
case. When we say that helicity amplitudes are measured or predicted this always im- 
plies a determination up to an overall constant complex factor. The usual assumption 
that all the helicity amplitudes are relatively real (cf. and references quoted there) 
reduces the numbers of parameters to four (three) for the decays involving each D**(l + ) 
(D**(2 + )). Each eliminated phase leaves, however, a sign ambiguity. A further reduction 
of the number of parameters requires a dynamical theory or model. We present absolute 
predictions for the angular decay distributions (no free parameters) from Heavy Quark 
Effective Theory (HQET) described in 0, g, 

The plan of this paper is as follows. In Section 2 we write down general formulae for 
decay distributions following the classical method of Jacob and Wick [JIH] . We introduce a 
diagrammatic representation of the decay formulae in order to keep track of the numerous 
terms. In Section 3 we use statistical tensors (see []T5| and references quoted there) to 
simplify the formulae and their diagrammatic representation. The discussion in Section 
2 and Section 3 is quite general and applies to simple decays as well as to cascade decay 
chains of arbitrary length. In Section 4 we apply our algorithm to the decays in eq. 
([]]). As an illustration we show how these formulae give an absolute prediction for all 
angular distributions when supplemented by the predictions of HQET. In Section 5 we 
discuss some problems relevant to the practical analysis of data. Section 6 contains our 
conclusions. 



2 Diagrammatic representation of decay formulae 

Writing down the joint angular distributions for the decays in eq. ([I]) is in principle very 
simple. In practice, however, it is cumbersome to keep track of the six spin reference 
frames and of the many summations. Therefore, we propose a diagrammatic method 
which eliminates most of the effort. In the present section we introduce our method by 
applying it to the (extended) approach of Jacob and Wick ]T5|. 

Let us generalize the decay process in eq. (JTJ) as follows. A parent particle (e.g. 
B) decays into two first generation particles (e.g. B — > D** + W). Each first generation 
particle may, but does not have to, decay into two second generation particles (e.g. D** — > 
D*tt and W — > Iv). Each second generation particle may, but does not have to, decay 
into two third generation particles (e.g. D* — > Dir). The process may be continued. 

Each of the iV two-body decays (N = 4 for process in eq. flU)) is described in the 
rest frame of its decaying particle. The masses of all the particles are considered to 
be known, therefore, the momenta of the decay products of each two-body decay are 
completely specified by the spherical angles 9, <fi of the momentum of the first decay 
product. The problem is to find the joint decay angular distribution W(9i, <pi, . . . , 9^, 4>n)- 
This distribution can be calculated as the trace of the joint spin density matrix of the 
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Figure 1: Building blocks for the diagrams 
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N + 1 final particles: 

w(e 1 ,...,<p N )= £ /'aI::!;:: ( 2 ) 

Ai,...,Ajv+i 

Let us denote the spin density matrix for a single particle by p\,. The indices A, A' 
usually denote helicities. Sometimes, however, helicity is not sufficient to define the spin 
state of a particle. An important example is the following: a virtual VT-boson with helicity 
zero can be either a spin one (vector) particle or a spin zero (time like) particle. In such 
cases it is understood that each index contains besides the value of helicity also the value 
of spin. This, of course, requires a suitable notation. For the virtual W boson one uses 
the subscript which is the helicity A = ±1, for the vector particle and t for the time like 
particle. In general, we interpret each index as a pair: (s, A) where s denotes the spin 
and A the helicity of the particle. If the particle spin is well defined we will omit the spin 
specification. For further use we ascribe to each element p\, of the spin density matrix of 
the parent particle the diagram shown in Fig |l|a. 

As a result of each two-body decay, the joint spin density matrix of the system acquires 
two more indices. Let us take the decay A — > B + C for definiteness . The description 
of this process is particularly simple when the spin reference frames shown in Fig. [| are 
used. For this choice of frames the axes za arid zb are parallel to each other and to 
the momentum of particle B, as seen from the rest frame of particle A . The axis zc is 
antiparallel to the axes Za, %b and parallel to the momentum of particle C as seen from 
the rest frame of particle A. The choice of the other axes is described below. 

Further we shall call the frames like the one shown in Fig [| for particle A decay 
frames. The decay frame for a particle is defined using the momenta of next generation 
particles. Frames like the ones shown in Fig. |2] for particles B and C will be called 
production frames. They are defined using the momenta of particles from the preceding 
generation. Each particle has a production frame, though the production frame may be 
related to a collision process rather than to a decay for the parent particle. Each particle 
that decays in a given process has a decay frame as well. The rotation which converts the 
production frame of particle % into the decay frame of the same particle will be denoted Ri. 
This rotation can be parameterized by the three Euler angles <f), 8, ip. According to our 
construction 9 and <p are the spherical angles of the momentum of the first decay product 
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Figure 2: Spin reference frames for the decay A — > B + C 



of particle i in the production frame of particle i. We choose the x axis in the decay frame 
of particle i so that the third Euler angle ip is zero. Geometrically, this means that for 
particle i the z production -a.xis, the z^ecay-axis and the Xrf ecaj/ -axis are all in one plane. The 
orientation of the x,i eca2/ -axis is such that going in the plane from the positive z pro d U ction- 
semiaxis past the positive Zd eca2/ -semiaxis one encounters first the positive Xd ecay -semiaxis. 
The ?/decay _ axis is chosen so as to make the x decay ,y 'decay > z 'decay frame right-handed. The 
l/prod«rfion _ axes for the decay products of particle % are chosen parallel to the 2/d eca2/ -axis of 
particle i as shown in Fig. ^|. The x pro d uction -axes for the decay products are chosen so 
as to make the corresponding production frames right-handed. Note that this makes the 
^production axis for particle B antiparallel to the x pro duction axis for particle C. 

The joint spin density matrix for particles B and C originating from the two-body 
decay of particle A is 

„sb,Ab,sc,Ac \ ^ saAa tj tj* fo\ 

Ps' B ,X' B ,s' c ,X' c ~ Ps / A ,X' A n SA,XA\s B M,so,Xo rl s' A ,X' A ;s' B ,X' B ,s' ,X' \°) 

saAa,s' a ,X' a 

The respective reference frames are defined in Fig. |2|, i.e. the spin density matrix 
for A is in the decay frame of A and the spin density matrices for B and C are in the 
production frames for B and C. H Sy \-_ are the helicity amplitudes where sa, Aa, sb, 
and 8c, Ac are the spins and helicities of A, B and C, respectively. The conservation 
of the ^-component of angular momentum implies that the helicities of the particles are 
related as follows 

Aa = A B - Ac and \' A = X' B - X' c (4) 

Relation (^|) corresponds to the diagram shown in Fig |3|a, where we have used the 
building block introduced in Fig |l]c. Summations over indices (helicities and spins) cor- 
responding to the lines without free ends are understood. In practice the spin density 
matrix of particle A is usually known in the production frame. Then in order to use 
formula (|3|) one must first rotate the spin density matrix of particle A according to the 
rule 

pfA = E pA'DxL (Ra)D^ x>a (R A ) (5) 

A, A' 
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Figure 3: Diagrams for the decay A — > B + C. 

a) Starting from the decay frame of particle A. 

b) Starting from the production frame of particle A. 



A^B + C 
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c>Ac 
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(b) 



: 'b> a b 



C' A C 



Figure 4: Diagram for the angular distribution of particle B from the decay A — > B + C . 




fax 



A^B + C 



where p and p denote the spin density matrices of particle A in the production frame and 
in the decay frame respectively. The diagram for the calculation of the joint spin density 
matrix of particles B and C in their production frames is shown in Fig ^b where one 
starts from the spin density matrix of particle A in its production frame. It includes the 
new building block shown in Fig |l|b which corresponds to the rotation of the spin density 
matrix of particle A from the production frame to the decay frame. 

As an example let us write down the angular distribution for particle B starting from 
the spin density matrix of particle A in its production frame. The diagram shown in Fig. 
|] corresponds to the formula 



w(9A) = Y,p°A> D x,*AM o)D( r>XA (4>,e,o) * 



H s ,x A ;X B ,x c H s , 



A/t ;An,Ar 



(6) 



where we sum over all spin and helicity indices and assume that the final decay particles 
have well defined spins. Here is defined by relations (^). We shall see in the next 
section how this formula can be simplified. 

3 Statistical tensors and simplified diagrams 

We propose to replace the elements of the spin density matrices used in the preceding 
section by the components of the statistical tensors 



'2J 



2s + 



-J2p s ^<s',X';J,-M\s,X> 

1 A, A' 



(7) 
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Figure 5: Simplified building blocks for diagrams. 




1 M 



tpJaJbJc 
Ma Mb Mc 



(cf. [15) and references contained there). Here < . . . | > are Clebsch-Gordan coefficients 
and the summation is over the values of A, A' corresponding to non vanishing Clebsch- 
Gordan coefficients i.e. — s < A < s, —s' < A' < s' and A' — A = M . With this 
replacement the building blocks of the diagrams are as shown in Fig. |||. Note that now 
each rotation introduces only one D-function, because formula fl5|) is replaced by (cf. JTi 
and references contained there) 



T M (s, s') 



E T^(s,s f )D J M , M (R) 

M' = -J 



The coefficients F^^/ 2 3 M from Fig [|b are calculated according to the formula 



JTlJl J2 J3 

r M 1 M 2 M 3 



(2J 1 + 1)(2J 2 + 1)(2J 3 + 1) 



E 



\\(2s A + l)(2s B + l)(2s c +l) x B x' g \ c x' c 

< s' A , \' A ; Ji, -Mi\sa, Aa >< s' B , \' B ; J 2 , -M 2 \s B , Xb > * 

< s c-> ^ci ^3) — M 3 \sc, Ac > H SA! x A;SBt x BiSC> \ c H*, 



1' V -c' V e' V 

*A' A A' S S' A B' A C' A C 



(9) 



Here A^ and A^ are given by the constraints (|J). The dependence of the coefficients 
FmiM 2 3 m 3 011 the decay and on the spins s A , ■ ■ ■ , s' c has not been written out. Formula (Q) 
implies 



Mi = M 2 - M 3 



and the identity 



(10) 



M1M2M3 



S A ,s' A ,S B , s' B , S C ,S C ) = (-1) 



^-Mi-M2-M 3 ( S A> S B> S c , Sc) 



which guarantees that the decay angular distribution is real as it must be. Taking the 
trace over a pair of indices of the joint spin density matrix is equivalent to coupling the 
two helicity indices to J = according to (|7|), then multiplying the result by \/2s + 1 
and finally summing over all the possible spins s of the outgoing particle. Denoting the 
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Figure 6: Two usefull identities in the diagrammatic form. 



J,M=0 




\J=0 



2J ^F{J-s,s>) 



(b) 



Figure 7: Simplified diagram for the decay A — > B + C. 




summation over a pair of indices by a corner without an outgoing line, diagrammatically 
one has the identity as shown in Fig. |6|a. 

As a simple example let us reconsider the decay A — > B + C '. The corresponding 
diagram is shown in Fig. [7] and replaces the diagram shown in Fig. [|. The formula for 
the angular distribution now reads 



J.M,s,s' 



(12) 



where we have used the well-known identity relating the rotation functions to spherical 
harmonics: 



D J MQ {4>,e^) 



2J+ 1 



(13) 



and the definition of the coefficient 18 



F(J;s,s') = N{-1 
< s, A - \';J,0\s',\ — X' > H, 



An; 



2s + 



XX' 



(14) 



where N is a suitable normalizing factor. The general identity between a box with one 
ingoing line and the coefficient F(J) is shown in Fig|6|b. Let us add that the coefficients 
F(J) can also be defined and evaluated for more than two body decays (cf. [0, |L9| and 
references given there). Thus our analysis can be easily extended to cascades, where one 
or more of the final decays is into more than two particles. 

Let us note that a circle in a diagram yields a factor Y^* whenever the index M of the 
line to the right of the circle equals zero. Because of identity ([10]) this is always the case 
for the circles preceding the final decays. Formula ( |i~2f) is clearly simpler than formula 
(M). Moreover, it implies the relation 



^MjE^r^^j; s,s') = J w{e,4>)Y^{e,4>)dn 



(15) 
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Figure 8: Simplified diagram for the cascade in eq. (1). 




which shows that each number a(J, M) can be obtained from the experimental angular 
distribution W(6, 4>) and that the set of numbers a(J, M) contains all the information 
obtainable from the joint angular decay distribution. 

One could perhaps object that the statistical tensors are less familiar than the spin 
density matrices. However, only the statistical tensors for the parent particle occur explic- 
itly in our formulae. In practice, the parent particle has either spin zero or spin one-half. 
In the former case there is only the J = statistical tensor 

T ° = 1 (16) 

In the latter, the allowed values for J are J = 0, 1 and the necessary components of the 
statistical tensor are 



T ° = ± (17) 

TL = TP^ ±h (19) 

The relation of the components of the statistical tensor to the elements of the spin density 
matrices can be seen to be very simple. In the following section we apply our algorithm 
to the decays in eq. (jl]). 

4 Joint angular distributions for the decay cascade 

B -> L>**(^ D*{^ Dn) + ?r) + Iv) 



The diagrammatic representation of the joint decay distributions for the decays in eq. ([!]) 
is shown in Fig 18]. The corresponding decay formula reads 
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Re 



2 2s D »* J 12 J 13 

w(6 l ,<f> l ,e D .,<f> D .,e D ,<i> D )= E E E E E ' 

J 3 =0,2 J 2 =0 Ji=0 A/=0Afi=-Ji 3 

a(Ji, J 2 , J 3 , M, Mi)D^} Ml ((f) D * , 9 D * , 0) 

(2 - Wm i0 ) (20) 



where the identity QTTp has been used. denotes the smaller of the two numbers Jj 
and Jj. The prime in the last sum means that only nonnegative Mi are to be included 
for M — 0. The coefficients a( Ji, J 2 , J3, M, Mi) are given by the formula 



a(J u J 2 ,J 3 ,M,M 1 )=F D *(J 3 )F$Xo E F^(s w , s' w )F w (J 2 ; s w , s' w ) (21) 



=0,1 



Each of the coefficients a(Ji, J 2 , J 3 , M, Mi) is independently measurable because for- 
mula (|20| ) implies the relation 

9 7 _i_ 1 /■ 

o(Ji, . . . , Mi) = J dQ l dn D ^ D W(9 h . . . , <j> D ) * 

'Djf Ml (<l> D *,e D .,O)Y&{0i 1 </>i)Y*(0 D , (f> D )] (22) 

The integral is just the average of the expression in square brackets over the experimental 
joint decay distribution W(9i,...). In principle, the averaging should be performed at 
fixed PF-boson mass \fq I . In practice, one would probably use some model for the in- 
dependence of the helicity amplitudes in order to be able to use all the data simultaneously. 
The coefficients Fqj^^(sw, s' w ) correspond to the decay B D** + W and the coefficients 
Fmm.0 to tne decay D** — > D*n. The two possible spins of the W-boson s w , s' w = 0, 1 
have been taken into account. For each of the two decays the spin state of one particle 
(A = ±1,0, t for the W and A = ±1,0 for the D*) and angular momentum conservation 
fix the spin states of the other two particles participating in the corresponding two-body 
decay. Therefore, we shall use the notation H\ and h\ for the helicity amplitudes of the 
two decays, respectively. We normalize the angular distribution W(9i, . . .) so as to get as 
simple formulae as possible. In some of the following expressions overall factors will be 
dropped. Note also that because of the equality between the subscript of Y^{9i,4>i) an d 
the first subscript of Dj} M ((f) D *, 6 D », 0), the azimuthal angles <\>\ and 4>d* always occur in 
the combination <fii + 4>d* — ft ~ X where x is the angle between the half planes bounded 
by the zsdecay axis and containing the momenta of the lepton and of the D*, respectively. 
Thus the distribution W(9i, . . .) is effectively given in term of five variables instead of six. 



The coefficients Fi(J) have been reviewed and listed in ref. |15|]. The nonzero values 
needed for the present case are 



F D *(0) = (23) 
^*(2) = -\h£z (24) 
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*W(0;0,0) = \/^Y^ (25) 




F W (1;1,0) = -F w (l;0,l) = -y — — (26) 
MO; 1,1) = ^ (27) 

Mi;M) - (28) 



^W(2; 1, 1) = (29) 

where lepton mass effects (m; 7^ 0) are taken into account by introducing the ratio ( q 2 is 
the invariant mass squared of the lepton pair) 

The two possible decay cases W 7- — > l~ + v\ and P^ 4 " — > / + + ^ are both included in eqs. 
(H) where only 1, 1) in (gj) is affected (± for W /± ). It is seen that the terms with 

J3 = 1 vanish in the expansion (|20|). 

Let us consider the coefficients F^m^ for the decays _D**(1 + ) — > D*7r. Parity conser- 
vation implies /ii = Therefore, the dynamics of the decay D**(l + ) — > D*7r introduces 
two free parameters (the overall normalization is not counted) into the joint angular dis- 
tribution, which may be combined into one complex parameter r], if we put 

h = V h i (31) 

According to HQET rj = 1 for D**(l+, j = \) -> D*tt and r] = -2 for j = 

|) -> D*ti. From @ one has (see e.g. [13,14]) 



(32) 

« = FSSS = #(1-H 2 ) (33) 

(34) 

(35) 

(36) 
(37) 

where the overall factor \hi\ 2 has been dropped. It can be seen that (Ji, M±) = (1,0) is 
forbidden. 

For the decay D**(2 + ) — > D*n parity conservation implies h\ = — and ho = 0. For 
the coefficients F^^ lQ one finds 

« - 7^ (38) 



pooo 

^000 ~~ 


^(2 + H 2 ) 


p020 

^000 ~~ 


— 1 — \v\ 

3 v 1 /1 


F 120 _ 
^110 ~~ 


Im(rj) 


F 220 

^000 — 


la+z\v\ 2 ) 


p220 

^110 ~~ 
P 220 

r 220 ~~ 


Re{rj) 
1 
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« = (39) 



21 



« = ^ (40) 
« = (41) 



21 

/3 



^ 2 2 ° ±2 o = -\l (42) 



« - WiS5 (43) 



« = -4a/4 (44) 



105 



^° ±2 o = ^ (45) 

where again the overall factor \hi\ 2 has been dropped. Thus the decay D**(2 + ) — > £)* + 7r 
introduces no free parameters into the joint angular distribution. One also notices that 
there are no contributions with J\ = 1, 3 or M\ = ±1. 

Finally, we must consider the coefficients Fq^^(sw, s' w ). Let us begin with the decay 
B -> D**(l+) + W. For s w = s' w = 1 we find 

«° = ^l#-i| 2 + l#o| 2 + |tfi| 2 ) (46) 
« = FZ = - l^-il 2 ) (47) 



^o 1 = ^(l^i| 2 + l^i| 2 ) (48) 



^on = \{H^Hl + H Ht) (49) 

« = « = ^(|i/- 1 | 2 -2|iJ | 2 + |iJ 1 | 2 ) (50) 

« = « = ^(|tfi| 2 -|tf-i| 2 ) (51) 

^on = ^oTi = li-H^ + HoH*,) (52) 

« = i(|^-i| 2 + 4|F | 2 + |^r) (53) 

*on = \(H^H* + H H*) (54) 

F °|| = tf^tf* (55) 



while for s w = Sw = 0: 



« = ^\Ht\ 2 (56) 
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« = -]]lm 2 (57) 

For the interference terms F^ M (0, 1) = — F 0J ] l J_ M (l, 0). For s w = 1, s' w = 0: 



^000 — 




(58) 


pOll 

^Oilil — 




(59) 


F 021 _ 
^000 — 




(60) 


F 021 

^Oilil — 




(61) 


5 into formula (p0|) one obtains the distribution 1^(0;, 


■ ■■A 



as a sum of 43 terms. Each of the 43 coefficients a(J±, J2, J^,M,Mi) is experimentally 



measurable according to formula (|T5|). These coefficients, however, are all functions of 
the same helicity amplitudes Hi, Ho, H t , hi, ho- Since the overall normalization and 
the overall phase of the amplitudes for each decay drop out from the formula for the joint 
angular distribution the number of real dynamical parameters that are to be extracted 
from the data is eight. The problem can be seen to be grossly overdetermined. Never- 
theless, each measured coefficient a(Ji, Ji, J^,M,Mi) can be used to reduce the error of 
measuring the helicity amplitudes. 

For the decay B -> D**{2 + ) + W and s w = s' w = 1: 

F Z = -^(l^r + l^r + l^l 2 ) (62) 



iffi = i(l^| 2 -|^| 2 ) (63) 



30 

F™ = -±=m\ 2 ~ \H-i\ 2 ) (64) 
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« = -4^(l^| 2 + 2|^or + |^-i| 2 ) (65) 



42 

^ooo = 4^(l^i| 2 - 2 l^o| 2 + |^-i| 2 ) (66) 



30 

*S8 = -^(l#ii 2 -l#-ii 2 ) (67) 

*E? = ^=(H H* - H.iH*) (68) 

« = -^={\H^-A\Ho\ 2 + \H^) (69) 

« = fj^iHl (71) 
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« = -^(2m 2 -3\H \ 2 + 2\H^) 
V35 



/ 3 

17041 . / / TT TT* TT TT*\ 

-^011 - _ Yii^ 1_ 0; 



*%2 = -2^(|^| 2 + 3|ifo| 2 + |^i| 2 ) 


p042 _ - TT TT* 

*022 - _ 7/7 1 1 



while for svk = s'm = one has 



pOOO _ / x I rr 1 2 
^ooo - Vcl^tl 



/2 

17020 / I rr 1 2 

-^ooo - -y^l-nd 
F 4 = 3,/I|// t ' 2 



35 

Finally, for the interference terms F^ M (0, 1) = — Fq£]^_ m (1, 0) and for sw = 1, s' w 



p001 _ Arj TT* 

^ooo - y s H ° H t 

77021 _ J u u* 

-^000 - -yj^o^t 

p 021 — J ^ u n* 

t Q±i±i - — V 14 Tl * 

6 

^000 - T/fo 



p041 /3 o , , 

*0±1±1 - V Tj H Tl H t 



72) 
73) 

74) 

75) 

76) 
77) 



78) 

79) 

80) 
= 0: 

81) 

82) 

83) 
84) 

85) 



In this case the expansion (|20| ) consists of 48 terms. The 48 coefficients a(J\, J2, J3, M, Mi) 
depend on six real parameters corresponding to the helicity amplitudes H 1 , H_i, H ,H t . 

In order to reduce the number of free parameters in the joint angular distribution 
to below eight in the D**(l + ) case and to below six in the D**{2 + ) case, it is necessary 
to know something about the dynamics of the decay processes involved. The HQET 
predictions for B -> D**(J P ) + W (cf. @, [TJ, 0, [TJ) are given by 



H t (l+) ~ H (2 + ) ~ J^±1(M B -M D ) (86) 



13 



H (l + ) ~ #*(2+)~ s j ( ±Jl( K M B + M D ) (87) 
1. 



#±i(l + ,j = -) 2if ±1 (l + ,j = -)~±vWl-v / ^T (88) 

#±i(2 + ) ~ ^-(v^TTtv^I) (89) 

where u is the velocity transfer variable uj = V\V^ as usual. 

In these formulae the symbol ~ means equal up to a factor independent of A. Factors 
dependent on J p and j have been dropped. Thus, these formulae should not be used 
to compare different decays. For given J p and j, however, they yield all the amplitudes 
H\(J p ,j) up to a constant factor and consequently are sufficient to calculate all the 
coefficients F^jlf from the formulae given in this section . Combined with the values 
of the parameter rj given previously, they give absolute and complete predictions for the 
joint angular distributions of all the cascade decays in eq. ([[]). 



5 Applications 

Once one has measured the coefficients a(Ji, J 2 , J3, M, Mi) it is possible to determine 
(up to a common complex factor) all the helicity amplitudes H\ and h\. In fact the 
problem is grossly overdetermined, so that many strategies and many consistency checks 
are possible. In practice the choice will be dictated by the available data. For the D**{2 + ) 
it is necessary to know at least six ratios of the coefficients a(Ji, J 2 , J 3 , M, Mi) in order to 
find all the amplitudes H\. This reduces to four ratios when the lepton mass is neglected 
and to three ratios, if the amplitudes are assumed to be relatively real (sign ambiguities 
remain). Neglecting both the lepton mass and the relative phases, one is left with two 
real numbers to be determined. In the D**(l + ) case one may be interested either in the 
amplitudes H x or in the amplitudes h\. We show below how the two problems can be 
separated from one another. 

Let us consider the ratios of the coefficients a(Ji, J 2 , J3, M, Mi) with common values 
of Ji, J 2 and M. In such ratios the amplitudes H\ cancel. If the lepton-side information 
is not used, i.e. taking J 2 = M = 0, we find for the D**(l + ) cases: 



a(0,0,2,0,0) 
a(0,0,0,0,0) 



2 


7/| 


2 - 1 


V5 


7/| 


2 + 2 



(90) 



a(2,0,0,0,0) = yg 

a(2,0,2,0,2) 3 Vl " ' K ' 

a(2,0,2,0,0) l H , o| , 2 



a(2,0,2,0,2) 3 
0(2,0,2,0,1 



o(2,0,2,0,2) 



o( 1 + 2 I^D ( 92 ) 
Re{rj) (93) 



These ratios yield the complex parameter 77 up to the sign of Im(rj). Since one expects 
theoretically that Im(rj) is close to zero, the determination of this sign may be difficult. In 
principle it can be deduced from the sign of the coefficient a(l, 1, 2, 0, 1). The vanishing 
of Im(rj), on the other hand, is easy to test, because it implies the vanishing of five 
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coefficients. Integrating the decay distribution over both dVti and dQo, i.e. considering 
only the decays D**{\ + ) — > D* + ir, one is left with two coefficients: a(0, 0,0,0, 0) and 
a(2, 0, 0, 0, 0). Their ratio depends both on the amplitudes H\ and on the amplitudes h\ 
and, of course, is not sufficient to find either set. This is the problem pointed out for 
the amplitudes h x in ref. |7[|. Our solution is to look at the joint decay distribution of 
D** — > D*{-^ D + 7r) + 7r. Contrary to ref. [[/], we make no assumptions about the spin 
density matrix of the D**. The HQET predictions are obtained by substituting rj = 1 
and rj = — 2 in the case j = \ and j = §, respectively. 

In order to find the amplitudes H\ for the D**(l + ) case, we propose to look at the 
decay distributions integrated over dQc i-e. at the coefficients with J 3 = Mi = 0, and 
to compare pairs of coefficients with the same values of J\. Then the dependence on the 
amplitudes h\ drops out. We find seven equations for the six unknown parameters: 



a(0, 


1 


0,0,0) 


a(0, 





0,0,0) 


a(0, 


2 


0,0,0) 


a(0, 





0,0,0) 


a(2, 


2 


0,1,0) 


a(2, 


2 


0,0,0) 


a(2, 


2 


0,2,0) 


a(2, 


2 


0,0,0) 


a(2, 





0,0,0) 


a(2, 


2 


0,0,0) 


a(2, 


1 


0,0,0) 


a(2, 


2 


0,0,0) 


a(2, 


1 


0,1,0) 


a(2, 


2 


0,0,0) 



±(|#! 


12 _ 




H 


-1 


2 ) 


- AeRe(H H;) 


2(1 + 6) | 


#1 


2 


+ 




-1 


2 + 


Hot 


2 



1 1 - 2e 


#i 


2 -2|tf 


2 + 




2 


2v/5 1 + e 




2 + 1 




2 + 




2 



|# 1 | 2 + 4| r f | 2 + (if.xl 2 



6 



2^5 



4|/f | 2 + (Jf^l 



1 + e 




2 — 2\H ot 


2 + 




2 


1 - 2e 


#1 


2 + 4|# 


2 + 




2 



\/l5 + (]#! 


2 _ 




H_ 


i 


2 ) 


+ 8eRe(H Ht) 


1 -2e | 


Hi 


2 + 4|# 


2 + 




2 



3^ T(#-i#o - #o#i*) + 2e(H^H; + H t H^) 



1 -2e 

where we have used the abbreviation 

|-^ot| 2 : 



|# 1 | 2 + 4| r f | 2 + 



(94) 
(95) 
(96) 
(97) 
(98) 
(99) 
(100) 



3e 



1 + e 



(101) 



From these relations one can find all the amplitudes H\ (up to a constant factor) and 
make one consistency checks. The HQET predictions can be obtained by substituting the 
amplitudes H\(l + ) from formulae (|86D - flS8f). 

The analysis for the D**(2 + ) case is much simpler, because there is no problem with 
the amplitudes h\. The integration over dQi and dVLp leaves, however, also in this case 
too few coefficients (a(0, 0, 0, 0, 0), a(2, 0, 0, 0, 0) and a(4, 0, 0, 0, 0)) and the study of some 
joint decay distributions is necessary. From the three coefficients one can only find the 
ratio 



R 





\Hot 


2 




2 + | 




2 



(102) 



and make one consistency check. 



15 



In order to experimentally determine the amplitudes h\ one can also use other decays. 
For the decay B — > /)**(—> £)*(—>• Dir) + 7r) + 7r the joint angular distribution can 
be formally obtained from the formulae given in Section 4 by putting H t = 1, if±i = 
-f/o = 0,.FV(0;0,0) = 1 and omitting the angular Oi, ^-dependence. The number of 
parameters a(J ± , J 2 , J 3 , M, Mi) left is six for D**{\ + ) and eight for D**{2 + ). For the 
D**(2 + ) case angular momentum and parity conservation (for the D** decay) are sufficient 
to make an absolute prediction. For the D**(l + ) case the six parameters depend on 
the complex parameter 77. The same argument applies when the first generation pion 
is replaced by any other spin zero meson. Another example is the nonleptonic decay 
B — > D**(^ D*(—>- Dtt) + 7r) + p(— > 7T7r). The necessary replacements in the formulae 
of Section 5 are the following: (i) H t = 0, (ii) omission of the arguments sw, s' w , (iii) 
replacement of 61, 4>i by 9 n , <fi n corresponding to one of the decay products of the p and 
(iv) insertion of F p (0) = F p (l) = and F p {2) = ^ for F w (0; 1, 1), F w (l; 1, 1) 

and Fy/(2] 1,1). The problem is similar to the problem for the decays in eq. ([!]) when the 
lepton mass is neglected. The number of measurable parameters, however, is reduced to 
25 for the D**(l + ) and to 32 for the D**(2 + ). The p-meson can be replaced by any spin 
one meson decaying into two spin zero particles. 



6 Conclusions 

The joint angular distribution for each of the decays in eq. (|l|) depends on more than 
40 independently measurable parameters. The coefficients a(J\, J2, J3, M, M\) are a con- 
venient complete set and can also be used to express other parameters. For instance, 
the lepton energy Ei in the B rest frame is a linear function of cos^. Consequently, the 
average < E\ > can be calculated as a linear function of the coefficient a(0, 1, 0, 0, 0). The 
average < Ef > depends, moreover, linearly on the coefficient a(0, 2, 0, 0, 0), while all the 
averages < E'j 1 > for n > 2 can be expressed as linear functions of these two. The decays 
of D**-s into longitudinally polarized D*-s studied in correspond to a determination 
of the coefficients a(Ji, 0, J 3 , 0, 0). 

The coefficients a( J±, J2, J3, M, Mi), though independently measurable, are strongly 
correlated as a result of angular momentum conservation. Thus, when testing models 
by comparison with data it is preferable to compare with the predictions for the helicity 
amplitudes and not for the coefficients a( J\, J2, J3, M, Mi). Otherwise it is easy to misin- 
terpret implications of angular momentum conservation as successful predictions of one's 
model. In Section 4 we have given the predictions of HQET, which would be very inter- 
esting to check. It is also important to test the relative reality of the helicity amplitudes. 
The helicity amplitudes can, of course, be replaced by form factors (cf e.g. fl5fl ). 



From experiment one would like to know in each case the coefficient a(0, 0, 0, 0, 0), 
(which defines the overall normalization), and as many other coefficients as possible, 
because each new measured coefficient reduces the experimental uncertainties. It is im- 
portant to realize that single-step decay distributions contain little information. Thus it 
is crucial to look at the joint decay distributions of at least pairs of particles. As discussed 
in Section 5, the joint two-particle decay distributions are sufficient to find all the helicity 
amplitudes. Nevertheless, the study of the full triple decay distributions can significantly 
reduce the errors in the determination of the helicity amplitudes. 
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